DC Josephson Effect with Fermi gases in the Bose-Einstein regime 
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We show that the DC Josephson effect with uitracoid fermionic gases in the BEC regime of 
composite moiecuies can be described by a nonlinear Schrodinger equation (NLSE). By comparing 
our results with Bogoliubov-de Gennes calculations [Phys. Rev. Lett. 99, 040401 (2007)] we find 
that our superfluid NLSE, which generalizes the Gross-Pitaevskii equation taking into account the 
correct equation of state, is reliable in the BEC regime of the BCS-BEC crossover up to the limit of 
very large (positive) scattering length. We also predict that the Josephson current displays relevant 
beyond mean-field effects. 

PACS numbers: 03.75.Lm, 03.75.Ss, 05.30.Jp, 74.50.+r 

I. INTRODUCTION 

In the last few years several experimental groups have observed, close to a Fano-Feshbach resonance [l[ , the crossover 
from the Bardeen-Cooper-Schrieffer (BCS) state of Cooper pairs to the Bose-Einstein condensateJBECLof molecular 
dimers in ultra-cold two-hyperfine-components Fermi vapors of 40 K atoms 0,0,0] and 6 Li atoms [a, H BUI • Few years 
ago the AC Josephson effect 0, in atomic BECs was predicted [ll| and observed fl2|. AC Josephson oscillations 
in superfluid atomic Fermi gases have been considered theoretically by several authors jl3l. [lil [TBL [ill [TtJ ■ Recently, 
Spuntarelli, Pieri and Strinati [l|| have studied the DC Josephson effect [H, [l(| across the BCS-BEC crossover in 
neutral fermions by using the extended BCS equations: they have computed the current-phase relation throughout 
the BCS-BEC crossover at zero temperature for a two-spin component Fermi gas in the presence of a barrier by 
solving the coupled Bogoliubov-de Gennes equations (BdG) [lsj . 

In this paper we show that a simple nonlinear Schrodinger equation (NLSE) [l6|, [l^, H(J, HH is able to reproduce 
the Josephson results of Spuntarelli, Pieri and Strinati 18], in the BEC side of the BCS-BEC crossover, i.e. from the 
deep BEC regime up to very larg e (positive) values of the scattering length. This NLSE is equivalent to the equations 
of superfluid hydrodynamics [22| with the inclusion of a gradient term [TH, H3, HH • We demonstrate, in particular, 
that the gradient term is essential to obtain the correct current-phase Josephson relation. 

As discussed in ref (HI, the DC Josephson currents found by the Gross-Pitaevskii (GP) and the BdG formalisms 
are the same in the deep BEC regime, while for relatively large values of the scattering length the BdG equations 
lead to important deviations from the GP predictions. Here we will show that our NLSE formalism, equivalent to GP 
in the deep BEC regime, gives results which are in remarkably good agreement with the BdG method also for large, 
positive, values of the scattering length. 

We also find that the breakdown of superfluidity, which corresponds to the maximum Josephson current across the 
barrier, strongly depends on the bulk equation of state embodied in the superfluid NLSE. In particular, on the basis 
of the Monte ' Carlo equation of state [20j that includes beyond mean-field effects, we predict that the critical currents 
are smaller than those calculated so far [18| using mean- field theories. 

II. NLSE FOR SUPERFLUID FERMIONS 

Inspired by the density functional theory of Helium 4 [23[ and by the low-energy effective field theory of the Fermi 
gas in the BCS-BEC crossover 0, H3 , we have recently introduced [HI, [2(| HlJ a complex order parameter 
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to describe boson-like Cooper pairs of a two-component fermionic superfluid made of atoms of mass m in the BCS-BEC 
crossover (22|, where n(r,t) is the local fermion density and 0(r,t) the local phase. Here n(r, t) = rt|(r,£) + nj(r,t), 
with n-\ (r, t) = n|(r, t). The normalization of ^(r, t) is such that 
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N being the total number of fermionic atoms. Notice that although ^(r, t) should not be interpreted as the condensate 
wave function [2y|, 8(r,t) represents the condensate phase 27], [28[. The local superfluid velocity v(r,t) = v-f(r,t) 
Vi(r,t) is thus related to the phase 9(r, t) by the equation 
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The nonlinear Schrodinger equation (NLSE) that satisfies Eq. ([3]) and reproduces the equations of superfluid hydro- 
dynamics [22j in the classical limit (K — ► 0) is given by 
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where U(r) is the external potential and (J,(n, ap) is the bulk chemical potential, i.e. the zero-temperature equation 
of state (EOS) of the uniform system, which depends on the fermion-fermion scattering length ap. In fact, by using 
Eqs. Q and ©, the NLSE can be written as 
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is a quantum pressure term containing explicitly Planck's constant %. This term can be viewed as a gradient correction 
in the density functional theory [23] or the next-to-leading correction in a low-energy effective field theory [25j . It is 
important to stress that in the deep BEC regime (ap — > + ) from Eq. ^ one recovers the familiar Gross-Pitaevskii 
equation for the Bose-condensed molecules made of paired fermions, where 
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with add(a F ) the dimer-dimer scattering length, which depends on the fermion-fermion scattering length ap. 
Over the full BCS-BEC crossover the bulk chemical potential can be written as [22j 



Kn,a F ) = ^(3, 2 n) 2/3 (f(y) - V -f(y) 



where f(y) is a dimensionless universal function of the inverse interaction parameter 
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with k F = (S^n) 1 / 3 the Fermi wavenumber and e F = h 2 k F /(2m) the Fermi energy. One can parametrize f(y) as 
follows: 



f(y) = u\ — ol-2 arctan a.z y 
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where the values of the parameters Oi, , «3 , (3\ , @2 , reported in Ref. [2(|, are fitting parameters based on asymptotics 
and fixed-node Monte-Carlo data [29|]. We call Monte-Carlo equation of state (MC EOS) the equation fj, — fi(n, ap) 
obtained from ([9]) and (fTTjl . Notice that Eq. ((4]) with Eq. (fTTj) has been recently used to describe density pro- 
files, collective oscillations and free expansion in the full BCS-BEC regime, finding a very good agreement with the 
experimental data [H, 0, [H, H3| • 

Within the mean-field extended BCS theory [28|, [3(| , the bulk chemical potential \i and the gap energy A of the 
uniform Fermi gas are found by solving the following extended Bogoliubov-de Gennes equations [26], |30( 
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From these two coupled equations one obtains the chemical potential /i as a function of n and of in the full BCS- 
BEC crossover (see for instance Ref. 26]). Note that, contrary to the MC EOS, this mean-field theory does not predict 
the correct BEC limit: the molecules have scattering length add — 2ap instead of the value add = 0.6af predicted 
by four-body and MC calculations [22T ]. We call mean-field equation of state (MF EOS) the equation a = /i(n, ap) 
obtained from Eqs. (fT2")l and (|13[) . Clearly this MF EOS is less reliable than the MC EOS, as shown in a recent study 

0- 

As previously stressed, the NLSE ([4| describes quit e accurately static properties and low-energy collective modes 
of oscillation in the full BCS-BEC crossover ■ fltil.ll9l.l20l]. but does not take into account the effect of pair breaking. 
In fact, Eq. ([!]) is reliable if the collective-mode wavelength A is such that A ^ £, where £ is the healing length of the 
superfluid. Recently, Combescot, Kagan and Stringari [31| have suggested that 

e = — (14) 

mVcr 

where v cr is the critical velocity of the Landau criterion for dissipation [13, [3l|. According to Combescot, Kagan 
and Stringari (3lj . in the BEC regime of bosonic dimers, and in particular for y > y c = 0.08, the critical velocity v cr 
coincides with the sound velocity, i.e. 
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Instead, for y < y c — 0.08, i.e. also at unitarity (y = 0) and in the BCS regime (y < 0), the critical velocity v cr is 
related to the breaking of Cooper pairs through the formula 
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where | A | is the energy gap of Cooper pairs [3l| . 

In the following section we shall show that the NLSE, Eq. ((4j, can be used to describe quantitatively the Josephson 
effect, but only on the BEC side of the BCS-BEC crossover, i.e. for y > y c = 0.08. For y < y c = 0.08 instead, pair 
breaking plays an essential role in the breakdown of the superfluid Josephson current [j| [3l[ . 



III. DIRECT CURRENT JOSEPHSON EFFECT 

We apply the NLSE (HI) to study the direct- current (DC) Josephson effect We consider a square- well barrier 

^( r ) I elsewhere 
which separates the superfluid into two regions, and assume a stationary solution 

*(r,t) - $(r) e l6 ^ e ~ a ^/ h (18) 
with constant and uniform number supercurrent 

J = n(r)v(r) = 2$(r) 2 — V0(r). (19) 
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From the previous equation it follows (V0) 2 = m 2 J 2 /(ft 2 $ 4 ) and also 

$(r) = 2/2 $(r) (20) 
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With the purpose of comparing our results with those obtained by Spuntarelli, Pieri and Strinati [l8j], we use the 
MF EOS of a{n{v),ap) and solve the above equation by imposing a constant and uniform density n far from the 
barrier region: 

Sto-yf for M-00 (21) 
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FIG. 1: Upper panel: scaled density profile n/n in the z direction (orthogonal to the barrier). The barrier width is L = 4/fcf. 
We consider three values of the energy barrier: Vo/tF = 0.025 (solid line); Vo/ef = 0.1 (dotted line); fo/eF = 0.4 (dashed line). 
Lower panel: phase difference A8 across the barrier for the same values of Vo/tF. 



We integrate Eq. (|20p on a ID mesh in real space, in an interval [0, z max ], using an imaginary time method, as 
described in [32j], and determine <&(]") by fixing the parameters of the barrier (which is located at z = z max /2), the 
uniform density n and the scattering length ap . To compute the current-phase relation we proceed as follows. The 
phase difference across the barrier can be obtained from Eq. (QI 
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We choose a value for A9 and compute J, from the above equation, at each iteration in imaginary time using the 
actual density profile $(z) 2 . The updated value of J is then inserted into Eq. (|20|) for the next iteration in imaginary 
time. We stop the calculations when convergence is achieved, i.e. the density profile does not change anymore between 
two consecutive iterations. 

In the upper panel of Fig. 1 we plot the scaled density profile n(z)/n of the fluid calculated for three different 
values of the energy barrier height. The figure shows that by increasing the energy barrier height Vo, the dip in the 
profile n(z) is enhanced. In the lower panel of the same figure we display the corresponding local phase difference 
A6(z) obtained from Eq. (fJU). These results are obtained by using the NLSE with the MF EOS. 

The calculated relationship between the current J and the phase difference AO is shown in Fig. 2 for two values 
of the interaction parameter y. Here our NLSE results (solid lines) are compared with the ones of Spuntarelli, Pieri 
and Strinati [HI , obtained by solving the full set of Bogoliubov-de Gennes equations for the quasi-particle amplitudes 
in the presence of the barrier (|17|) . For both values of y shown, the agreement with the BdG theory is remarkably 
good. For comparison, we also show in Fig. 2 the results obtained by solving the GP equation, to which our NLSE 
is expected to reduce in the deep BEC regime (y 1). Interestingly enough, we find that y = 3 is already in the 
regime well described by the GP formalism. At y = 1, however, large deviations of the GP curve with respect to the 
BdG (and NLSE as well) ones are found. 

In Fig. 3 we plot the maximum J™ ax of the current J as a function of the inverse interaction parameter y = 
l/{kpap), and compare our data (curves) with the results of Spuntarelli, Pieri and Strinati (symbols) [l8| . Figure 
3 shows that the NLSE reproduces the DC Josephson results of Ref. [HI in the BEC regime remarkably well, from 
the deep BEC regime (|/ > 1) up to very large positive values of the scattering length (y In the BCS regime 
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FIG. 2: DC Josephson current J vs phase difference A8 for two values of the inverse interaction parameter y — l/(kF(iF)- 
Barrier parameters: L — &/kF, Vo/cf = 0.1 (curves corresponding to y = 1), L = 5.3/fc.F, Vo/ef = 0.05 (curves corresponding 
to y = 3). Squares: Bogoliubov-de Gennes calculations of Ref. [19. [331] . Solid lines: NLSE with MF EOS. Dotted lines: GP 
results. 

(y < 0) the NLSE predictions are instead expected to be completely unreliable. This is hardly surprising since the 
superfiuid NLSE completely neglects the effect of pair breaking. 

Fig. 2 shows that deviations from the GP results of the calculated current-phase relation are found for the case 
y = 1. To further investigate to which degree the NLSE reduces to the GP case, in Fig. 4 we compare the predictions 
of the NLSE and those of the GP equation for the maximum of the Josephson current as a function of y. It appears 
that for relatively large values of the interaction parameter the NLSE results dramatically deviate from the GP results 
(while being in good agrement with the BdG calculations). In fact, on the basis of the Eq. (fl~5|) and the pair-breaking 
argument of Combescot, Kagan and Stringari [31], our NLSE should be accurate from y <C 1 (deep BEC limit) up to 
y = y c = 0.08, thus also for values of y very close to the unitarity limit (y = 0). 

Thus, the NLSE might represent a viable alternative, in the whole range of positive scattering lengths (and especially 
where the GP equation is not reliable, as shown in Fig. 4), to the much more computationally expensive BdG approach. 
This is especially true in the case of 3D geometries, where the BdG method could be prohibitively costly. 

As expected in our calculations we recover the Josephson equation 

J = J sin(A0) (23) 

in the regime of high barrier (weak-link). 

We observe that, contrary to our NLSE, the classical hydrodynamics equations of Fermi superffuids [22], i.e. Eqs. 
([5]) and ([6j) with Tqp — 0, cannot be used to study the Josephson effect. This is shown in Fig. 5 where we plot 
the current-phase diagram and compare the NLSE results (dashed lines) with the ones obtained using the classical 
hydrodynamic equations (dotted line). It clearly appears that the Josephson relation (|2"3"|) is violated (the dotted curve 
does not goes to zero at 7r) if we omit the gradient term. Moreover, the predicted current values are very different 
from the NLSE results. 

The results discussed up to now are based on the use of the MF EOS, since we were interested in assessing the 
reliability of the NLSE by comparing its results with those obtained by solving the Bogoliubov-de Gennes equations 
[l8l |. equivalent to a MF treatment. Of course, to give useful predictions to be compared with experiments one must 
use the MC EOS. In fact, relevant beyond-mean-field effects in the BEC side of the BCS-BEC crossover have been 
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FIG. 3: Maximum Josephson current J™ 11 * vs inverse interaction parameter y = l/^kpciF) in the BEC regiion (y > 0). Solid 
curve: J(J lax based on pair breaking in the BCS regime [Isl , |3l| . Other curves: superfluid NLSE. Symbols: Bogoliubov-de 
Gennes calculations of Ref . [18l |. Four values for the energy barrier height Vb/fF are considered: 0.025, 0.10, 0.2, 0.4 (from top 
to bottom in the figure). The width of the barrier is L = i/kp. 



predicted and observed for the density profiles and also for collective oscillations @, H3, [HJ . In a recent experiment of 
Miller et al. Q critical velocities have been observed in an ultracold superfluid Fermi gas throughout the BCS-BEC 
crossover. These critical velocities, determined from the abrupt onset of dissipation when the velocity of a moving 
one-dimensional lattice is varied [§[, are the analog of the maximum Josepson current J™ aa: . In Fig. 6 we plot J™ ax as 
a function of the inverse interaction parameter y — 1/ '(^pop). &nd compare the results obtained from the superfluid 
NLSE by using the MC EOS (solid lines) or the MF EOS (dotted lines). The figure shows that, for a given barrier, 
the maximum current predicted by the MC EOS is appreciably smaller than the MF one for all values of y. Moreover, 
as also observed in the experiment of Miller et al. Q , there is a pronounced peak of J™ ax at unitarity (y — 0) , which 
is absent in the MF curves. It is important to observe that, at the crossover, beyond-mean field effects exist not only 
in the bulk equation of state. In fact, we have recently shown that at unitarity (y = 0) the superfluid NLSE of Eq. 
(|4|) must be modified with the inclusion of an additional nonlinear term [34j . Nevertheless, this beyond- mean-field 
term goes to zero for a large number of atoms. 



IV. CONCLUSIONS 



We have introduced a NLSE equivalent to the hydrodynamic equations of Fermi superfluids plus a gradient cor- 
rection. Both hydrodynamics equations and superfluid NLSE are known to reliably reproduce static properties and 
low-energy collective dynamics. The advantage of using the NLSE is that one can take into account also surface and 
shape effects, and these can be relevant for a small number of particles. In addition, the gradient term is essential 
to obtain the correct Josephson equation, as demonstrated in the present work. We have shown that in the study of 
the DC Josephson effect our NLSE works quite well at the right side (BEC regime) of the BCS-BEC crossover also 
for very large values of the scattering length, where the familiar Gross-Pitaevskii equation is instead unreliable. In 
particular, our results suggest that the superfluid NLSE is accurate from y <C 1 (deep BEC limit) up to y = y c = 0.08, 
i.e. very close to the unitarity limit (y = 0). 

We thank Andrea Spuntarelli and Pierbiagio Pieri for making available their data [l8| and Nicola Manini for useful 
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FIG. 4: Maximum Josephson current Jo lax vs inverse interaction parameter y — l/(kF0,F) obtained with the NLSE. Comparison 
between BdG (dots) and the GP approximation (dotted lines) is made. Vo/cf = 0.1 and L — 4/kF- 
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FIG. 5: Current J vs phase difference A#: without the gradient term (dotted line) the Josephson-equation is violated and the 
result are very far from the NLSE ones (dashed line). 
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FIG. 6: Maximum Josephson current J™ ax vs inverse interaction parameter y — l/(kFO,F) obtained with the NLSE. Comparison 
between MC EOS (solid lines) and MF EOS (dotted lines). Two values of the energy barrier height are considered: Vo/cf = 
0.025 (upper curves) and Vo/tF = 0.2 (lower curves). 
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